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Traction drives for automotive applications have a rather checkered history. After a brief period 
of acceptance during the first quarter of this century they disappeared in rapid fashion. The main 
reason for their initial acceptance was the simplicity with which a variable speed output can be 
obtained from a nearly constant input speed. Their rapid demise was due to the fact that, while 
simple in concept, the average life expectancy was very low indeed. With the advent of a more 
demanding motoring public, this short life expectancy became unacceptable, and other solutions 
were found to the speed match problem. In recent years traction drives have had a resurgence because 
of the energy crisis and the fact that a lot of the initial problems have now been, partially, at least, 
solved. 

In simple terms the basic elements in a traction drive are two rollers, pressed into nominal 
contact and rolled about their respective axes. Power is transmitted in the form of a shear stress 
across the contact area. A fluid is present to prevent initial surface scuffing damage and to provide 
for some form of cooling. The rolling motion of the disks draws this fluid into the contact zone, and 
a thin layer of this fluid will separate the actual contact area. It is also in this region where the torque 
is transmitted from one roller to the next, and it should not surprise one that the performance of a 
traction drive depends to a large extent upon the rheological properties of the fluid. Close 
examination of the fluid as it passes through the contact gap reveals that it experiences a sudden 
pressure pulse from atmospheric to possibly several gigapascals in a period of 1 to 0.1 msec. The 
shear stress that is transmitted from one disk to the other (about 10 percent of the normal stress) 
passes through this layer of fluid “trapped” in the contact and causes a shear. This, in turn, will lead 
to heat generation, and from simple calculations, temperatures in the center of the film can easily 
reach several hundred degrees centigrade. 

To study the rheological properties of the fluid under these conditions precludes the use of most 
of the conventional instruments used for steady-state measurement. In fact the only suitable type of 
instrument is a disk machine where most of the conditions are the same or similar to those in traction 
drives. From the resulting traction tests, certain models are inferred, and it is in this area where there 
has been a lot of activity recently. 

To the designer of traction drives, the traction behavior of the fluid under the severe conditions 
is of utmost importance because of the direct influence that it has on the efficiency, size, and life of a 
given drive. 


Previous Traction Investigations 

As mentioned in the introduction, there has been a lot of activity in the area of traction research, 
both in the past and recently. Notable contributions have come from Clark, et al., Hewko, Smith, 
Smith, et al., Johnson and Cameron, Niemann, and more recently Johnson and Roberts and 
Johnson and Tevaarwerk (refs. 1 to 8). Some of these investigations were strictly experimental, and 
aimed at obtaining traction-drive design data, while others were aimed at understanding the traction 
phenomena so that rheological models could be formulated. This latter research is, of course, 
ultimately aimed at relating fluid molecular properties to traction properties. Research by Johnson 
and Tevaarwerk (ref. 8), Daniels (ref. 9), Hirst and Moore (ref. 10), and Alsaad, et al. (ref. 11) is 
directed specifically towards this purpose. The reader is referred to an excellent review by Johnson 
(ref. 12) for further aspects of this topic. 


•Work done under contract DEN 3-35. 
tAppIied Tribology, Ltd. 


285 


Many of the rheological models derived so far have been isothermal because of the degree of 
complexity that thermal analysis introduces. Thermal effects are important. However, a simple 
method is sought to include them in the analysis. 

Current understanding of traction has led to traction models that describe the fluid shear 
behavior in terms of an elastic and a dissipative element. An analysis of traction-drive performance 
using such a model was done by Tevaarwerk (ref. 13). It showed that under certain conditions the 
prediction technique by Magi (ref. 14) can be used. The work of reference 13 has recently been 
further expanded by developing a simple method to correct for thermal effects due to spin (ref. 15) 
and by performing an overall thermal traction study (ref. 16). 

As with all models, however, their usefulness is severely restricted if inadequate input traction 
data are available to the designer. This is especially so if high-traction fluids are used that have not 
been tested extensively for use under the wide range of conditions that exist in modern, highly 
advanced traction drives. 

Nomenclature 




First used 

>1/ 

regression constants 

eq. (11) 

a,b 

semi-Hertz contact size for x,y direction, m 

eq. (1) 

Ci 

regression constants for traction variables 

eq. (14) 

Fx 

contact force in the x-direction, N 

eq. (5) 

Fz 

normal force on the contact, N 

eq. (5) 

G 

average fluid shear modulus for contact. Pa 

eq. (1) 

h 

EHL film thickness of the fluid in contact, m 

eq. (1) 


dimensionless longitudinal slip variable 

eq. (4) 

h 

dimensionless side slip variable 

eq. (9) 

h 

dimensionless spin variable 

eq. (6) 

Ja 

dimensionless traction variable 

eq. (5) 

k 

contact aspect ratio, = b/a 

eq. (4) 

^1,2 

coefficients for rigid/plastic model 

eq. (8) 

m 

initial slope for zero spin traction curve 

eq. (1) 

P 

Hertzian contact pressure, GPa 

eq. (11) 

S 

dimensionless slip grouping 

eq. (9) 

U 

rolling speed of the disks, m/sec 

eq. (4) 

AU 

longitudinal slip velocity of the disks, m/sec 

eq. (4) 

AV 

side slip velocity of the disks, m/sec 

eq. (9) 

Z 

dimensionless spin grouping 

eq. (15) 

e 

speed pole parameter 

eq. (7) 

e 

temperature of the inlet fluid, °C 

eq. (11) 


peak traction coefficient for a given test 

eq. (2) 

Tc 

critical shear strength of fluid. Pa 

eq. (2) 

CO 

spin velocity on the contact, rad/sec 

eq. (6) 

Superscripts 

/ 

nondimensional regression variable 



denotes averaged variables 



Traction Predictions 

The principal traction information for the designer comes from the so-called traction curve. This 
curve is obtained by sliding and rolling two disks relative to each other and measuring the resulting 
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traction force. The data are normally plotted in dimensionless form as traction coefficient versus 
slide to roll ratio. Figure 1 shows such a typical curve. The exact shape and features of this curve 
depend on several external variables such as type of fluid, contact pressure, rolling speed, inlet 
temperature, and the kinematic conditions that exist in the contact. These kinematic conditions may 
consist of spin, side slip, and longitudinal slip, or a combination of these. It is, of course, not 
practical to experimentally measure the traction under every single variation of these conditions, so 
traction models are resorted to. With a suitable traction model and the knowledge of a few of the 
fluid rheological parameters, traction under varying kinematic conditions can then be predicted. 
Present-day traction models take into account the fluid and disk elasticity effect as well as the 
dissipative behavior of the fluid film. This model can be simplified for purposes of traction-drive 
analysis to a simple linear elastic spring in series with a Coulomb friction model (fig. 2). 

The method of operation of this model is as follows: Upon the application of a shear force, r, 
the linear spring element will extend. This extension is directly proportional to the stress at all levels 
of stress until a certain critical stress is reached. At this point the Coulomb friction is insufficient, 
and the extension in the spring element ceases. Sliding will now take place provided that the stress is 
maintained at t,.. Upon a reduction or reversal of stress, the elastic element will retract and recover 
the energy stored in it. The frictional work is, however, lost. (This frictional work will result in heat 
generation which, in turn, results in a local temperature rise. This temperature rise influences the 
level of stress, Tc, in a complicated manner but always reduces it.) The modulus, G, of the elastic 
element is related to the initial traction slope through contact geometry and the kinematic conditions. 

For simple slip traction curves, m and G are related as follows: 



h 

a^b 


( 1 ) 


where h is the average film^hickness, a and b are the semicontact dimension, and is the contact 
normal load. The quantity G is indicated here rather than G. This is to signify that we only deduce an 
average shear modulus. The reason for this is that we do not know exactly how the fluid shear 
modulus is influenced by the local pressure in the contact. (The extraction of G for the fluid is further 



Figure 1. - General features of a traction curve showing the initial 
traction slope m and the peak traction coefficient p. 


V\/V^- 


tc 


Figure 2. - Schematic of the elastic/plastic traction model, 
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complicated by the fact that in the traction slope a certain amount of elastic creep of the disk material 
is invariably present.) 

The frictional stress, is directly related to the peak traction coefficient, through the 
following equation: 




A 

irab 


( 2 ) 


Again, we are only extracting the average quantity because of the lack of knowledge about the 
variation of with pressure. (This limiting stress is temperature dependent. The temperature of the 
fluid is not constant but rises as more and more slip occurs on the contact. For simplicity, however, 
we often use the peak traction coefficient, if such a peak exists.) A detailed description of how this 
model can be used to predict traction under conditions of spin, side slip, and the usual longitudinal 
slip is given by Tevaarwerk and Johnson (ref. 17) and Tevaarwerk (ref. 13). The analysis described in 
these references is entirely isothermal in nature. 

To show how powerful this model is, consider the traction curves shown in figure 3 and reported 
by Tevaarwerk (ref. 18). These are curves obtained under a variety of conditions on two different 
traction testers and with two different fluids. Yet, by using an elastic/plastic traction model, all these 
curves can be reduced to fall on a single master traction curve (fig. 4). Each individual traction curve 
requires now only two parameters to describe it completely: the initial slope, m, and the peak traction 
coefficient, jji. This applies only to the simple slip traction curves, but the important parameters, m 
and /i, can be used to predict traction under any combination of spin and slip. 


Theoretical Traction Results 

The theoretical traction results, as calculated by Tevaarwerk and Johnson (ref. 17) and 
Tevaarwerk (ref. 13), for the elastic/plastic model are now presented. The reader is referred to the 
above two papers for further details. Although for a number of special cases analytical solutions 
exist, the general situation of combined slip, side slip, and spin was solved using digital techniques. 

Traction Due to Slip 

For simple dimensionless slip, J\, on the contact aspect ratio, k, the resulting traction, J^, is as 
follows: 


1 , 
* X 


2 tan-' S + 


2S 

1+S2 




where 5 = | — 

3 Vk 


( 3 ) 


and the dimensionless slip is 

, 3X m AU rr 
a H U 

and the dimensionless traction is 


( 4 ) 


■/ 4 =^ ( 5 ) 

This traction is a function of J\ and k, as shown in figure 5. It may be observed from figure 5 that 
contacts of low aspect ratio have a superior performance in terms of traction at a given slip value. 


Influence of Spin on Traction 

For the combined action of slip and spin the resulting traction is given in figure 6. Dimensionless 
spin is defined as 
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Fluid 

curve 

inlet 

Hertz 

speed 

aspect 

normal 

type 

number 

temp 

press . 

U 

ratio 

load 

(-) 

(-) 

(°C) 

(GPa) 

{ m/sec ) 

(-) 

(N) 

SANTO 50 

1 

32 

1.45 

20 

5 

5000 

" 

2 

69 

1.0 

30 

5 

1500 

'■ 

3 

29 

1.9 

20 

1 

1400 

*' 

4 

69 

1.0 

80 

1 

200 

TDF-88 

5 

30 

1.9 

10 

1 

416 


6 

70 

1.45 

30 

1 

185 

** 

7 

70 

1.0 

30 

5 

460 


8 

33 

1.45 

10 

5 

1400 


Figure 3. - Typical traction curves obtained on two fluids at a range of 
speeds, aspect ratios, inlet temperatures, and contact pressures. 


h = 


3tt 

T 


m 

y- 


coV^ 

u 


'iic 


( 6 ) 


No known analytical solution exists, and the results were obtained by digital integration. The 
influence of spin is to increase the slip, though this appears to be so only above a certain threshold 
value of spin. When the spin is sufficiently high, it is expected that elastic strain may be neglected and 
the rigid/plastic model may be used. 

The resulting traction due to the shear of a rigid/plastic-like material in the contact is shown in 
figure 7. The independent parameters are the aspect ratio k and the speed pole location e. The speed 
pole location parameter is identical to that used by Magi (ref. 14). It is related to the dimensionless 
parameters used in this analysis as follows: 
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Figure 6. - Dimensionless traction curves for the elastic/plastic traction 
model under combinations of spin, Jo, and slip, J-i , for an aspect ratio 
of /k = 1 . 



Figure 7. - Dimensionless traction curves for the rigid/plastic traction 
model at different speed pole locations. The speed pole location is given 
by E = J^/ J] • 
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(7) 


_ 7l _ AU 

Comparison with results obtained by Magi showed complete agreement. The results were 
obtained by numerical integration, and no known analytic expression exists for the traction in terms 
of the speed pole parameters. However, the curves shown in figure 7 do lend themselves to a simple 
empirical expression of the form: 

74=l-exp[-(eA:,)^2] (8) 

This expression can also be written explicitly for e. Coefficients k[ and ki are given in table I for 
various aspect ratios. The largest error between the theoretical and predicted traction, according to 
equation (8), occurs in the low traction region and is less than 7 percent of the peak traction 
coefficient. 

The traction curves and the influences of spin, as shown in figure 6, can be divided into three 
regions of influence. These regions are shown in figure 8 as a function of aspect ratio and 
dimensionless spin. The three regions represent: 

(1) Traction predicted by neglecting spin and by using the elastic/plastic model 

(2) Traction predicted by using the elastic/plastic model and by including spin 

(3) Traction predicted by using the rigid/plastic model and by including spin. 

It is in this latter region, where thermal effects become important, the methods outlined by 
Tevaarwerk (ref. 15) may have to be used to obtain more accurate traction predictions. 

Influence of Side Slip on Traction 

For the combined action of slip and side slip, the resulting traction for an elastic/plastic 
model are given by 


r M 

' 3 sVk 

where 

^ 3-ir m AV ^ 

■'2'T 71/ 

and 


_ \_ 
7T 


2 tan - 1 S + 


2S 

1 +52 


] 


(9) 


Figure 9 shows traction J 4 as a function of slip J\ at various levels of J 2 for an aspect ratio of 
k=l. The effect of side slip J 2 is to reduce the obtained traction in the working range. 

When the imposed side slip becomes sufficiently large, the elastic effects may be neglected, and 
the traction can then be obtained by the rigid/plastic model: 


/4 = 


Jj 


( 10 ) 


The reader might inquire as to the importance of being able to calculate the effect of side slip on 
traction. Side slip in traction drives results from misalinement of the axes of rotation of the principal 
shafts. Only very little misalinement is needed before 0.1 to 0.2 percent side slip is introduced. This 
can influence the longitudinal traction significantly, especially when little spin is present on the 
contact. 
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J3(apln) 


TABLE 1. -COEFFICIENTS 
FOR THE RIGID/PLASTIC 
TRACTION CURVE 
EXPRESSION 


Aspect ratio, 
k 

*1 

A:: 

0.5 

1.8 

0.9 

1.0 

1.9 

1.1 

2.0 

1.4 

1.5 

4.0 

1.1 

1.5 

8.0 

.8 

1.5 



k (aspect ratio) ^ b/a 


Figure 8. - Regions of influence of spin on the 
traction prediction of an elastic/plastic model. 
Dimensionless traction, is maintained at 0.75. 



J 1:0 10.0 100.0 


(Ji) SUP 

Figure 9. - Dimensionless traction curves predicted 
for combinations of slip, J-j, and side slip, J2» 
using the elastic/plastic model. 
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Fluid Traction Data 


Without the knowledge of the two fluid traction parameters, m and fi, the above models are not 
much use to the traction-drive designer. They do indicate traction trends when one or two parameters 
are varied, but no quantitative data can be obtained from them. To use the theory effectively one 
does need m and These are related to fundamental fluid properties as we have seen before, but how 
these fundamental properties tie in with molecular structure of the fluid is not yet fully understood 
and often one resorts to experimentally available data for the particular fluid of interest. 

To generate the experimental traction data for two fluids the Bearing, Gearing, and 
Transmission Section of NASA Lewis managed an extensive traction research under a Department of 
Energy program. This research was undertaken by Transmission Research, Inc. The main purpose of 
this research was to establish the traction parameters for two high-traction fluids at high speeds (10- 
to 80-m/sec surface velocity), at various contact aspect ratios (A: = 0.5 to 5), at inlet temperatures of 
30° to 70° C, and at contact pressures as encountered in traction drives. Also, the traction was 
measured under a variety of spin, side slip, and longitudinal slip combinations. 

The traction data obtained from this investigation is reported by Tevaarwerk (ref. 18). The two 
important parameters for the traction curves were extracted and correlated with the various 
independent parameters that were imposed. 


Traction Data Correlation 


Examination of the traction parameters obtained by the above investigation revealed that there 
are certain trends when one of the external parameters, such as pressure, speed, inlet temperature or 
aspect ratio, is varied. The trends are 

• Traction coefficient decreases with decreasing pressure. 

• Traction coefficient decreases with increasing speed. 

• Traction coefficient decreases with increasing temperature. 

• Traction coefficient decreases with increasing aspect ratio. 

Similar observations were found on the initial traction slope. By linearization of these trends the 
following equations may be written: 


/xoc/l j +AiP 


A^U 


-j 


(xocAi+A^k 


( 11 ) 


These equations can be used directly for correlation purposes; however, to reduce the influence 
of weighting, because of the physical magnitude of certain of the variables, nondimensional variables 
should be used. These can be obtained quite easily by using transformations of the type 


U' = 


V-L' 

1/ 


p-p e-e 

U = — 

p e 



( 12 ) 


where the prime indicates the transformed variable and the barred variable is the average magnitude. 
Because the trends are all equally likely for a given condition of speed, temperature, pressure, and 
aspect ratio, the following general equation may be expected: 


IJL = {Ai-i-A2P')iA-i + A4U')^^$+ -^^(A-j + Afjc') (13) 

Multiplication yields the following expression: 


M = Ci + C2/c' + C3P' + C4P'At'+ P+q|!-+C7^ +C^k'~+C^U' + C^oU'k' 


+ C„C/'P'+C,2f/'P'A:'+C,3^ +Cuk'^ TCigAr'P' ^ (14) 
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where the coefficient Cj to Cig are directly related to the coefficients Aj to Ag. 

Initial Traction Slope 

Equation (14) was also used on the initial traction slope data to obtain the coefficients C\ to C]6 
by a least squares regression. The resulting regression coefficients are shown in appendix A for the 
two fluids tested. For these regression coefficients all the traction slopes for a given fluid were 
nondimensionalized by using equation (12). Before one can use the regression results to predict the 
initial traction slope, the variables need to be transformed according to this equation. This requires 
the knowledge of the average quantities of the variables and these are given in appendix A. Thus 
equation (14) can also be used to predict initial slope in using the appropriate coefficients C\ to Cj^ 
given in appendix A. 

Traction Coefficients 

A similar analysis to the above could be done for the peak traction coefficient; however, it is also 
influenced by thermal effects due to spin. When spin is present, therefore, a further variation in the 
traction coefficient is possible. A similar treatment to the initial traction slope results gives the 
following expression for the traction coefficient: 


fi = C\ + C-^' + Cgk' + C^k' Z' + C^P' + Cf^' Z' + C-]P' k' + CgP' k' Z' + + C]q + Cn 

+ Cnk'y+C,g^ +C,^P'^ +C^^k’^ +C,(,k'P'y +C^,U’ +C^gU'Z' 

+ Ci<)U'k' +C2oU'k'Z' +C2\U'P' +C22U'P'Z' +C2iU'P'k' +C2AU'P'k'Z’ 


+ C25 ^ +C 26 U' ^ +C 2 jk’ ^ +C 2 sU'k’ ^ +C 29 P' y +C 2 oU'P’ ^ 
+ C2ik’P'^ +C22U'k'P'^ 


(15) 


where 


Z = 


u 


and 


Z 

The regression coefficients for the traction coefficient of the two fluids are given in appendix A. They 
can be used to predict the traction coefficient whether spin is present or not. 

It could be argued that a regression of the above type is not desirable because of the large 
number of coefficients involved. This is certainly true, and it is hoped that at some future time all the 
data can be reduced to two or three fundamental fluid parameters. However, such a model has not 
been sufficiently developed for inclusion here. 


Example Traction Predictions 

As an example of how the foregoing material may be used, let us take a look at some existing 
traction information and see how well the current methods will predict this traction. We will take 
some traction results as reported by Gaggermeier (ref. 19) and shown in figure 10. These data were 
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Figure 10. - Comparisons of measured and predicted traction under slip and 
side slip combinations. The following side slip was applied: ■ and curve 
A, 0 percent; a and curve B, 1.57 percent; • and curve C, 4.93 percent. 


obtained on Santotrac 50 at a peak pressure of 1.854 GPa, an inlet temperature of 50° C, and a 
rolling speed of 8.4 m/sec. The aspect ratio for the contact was 1.5. The three sets of data points 
shown are for different degrees of side slip. These data, then, are a combination of side slip and 
longitudinal slip. The traction prediction based upon the elastic/plastic model as presented here is 
given by equation (9). To predict the actual traction, we need to know the initial traction slope and 
the peak traction coefficient for the particular fluid under the experimental conditions. These two 
parameters may be predicted by using the correlated data as reported here. We need to use equations 
(14) or (15) together with the coefficients as listed in appendix A. For Santotrac 50 at the particular 
conditions, as applicable to the experimental data, we get 


m = 55.3 

fi=0.I06 


(16) 


These two values can then be used in predicting the traction from equation (9). The predicted traction 
is shown in figure 10 as solid curves. It may be observed that the predicted values are somewhat 
higher than the experimental data. The most probable reason for this difference is the influence of 
thermal effects, which tend to reduce the traction at higher slip values. To get an improved 
prediction, methods as outlined in references 15 and 16 should be used. 


Conclusion 

In this paper it has been demonstrated how the traction performance of a typical drive can be 
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predicted from the basic knowledge of two fluid parameters, the initial traction slope, and the peak 
traction coefficient. These two parameters are given for two high-traction fluids as a function of inlet 
temperature, contact pressure, surface velocity, aspect ratio, and, in the case of the peak traction 
coefficient, the amount of spin on the contact. The traction behavior of the two fluids is very similar, 
with only minor differences. Thermal effects are not treated in this paper but can play a major part in 
the traction behavior of a contact. Current research indicates that a fairly simple thermal correction 
gives very good results. Also, very recent research by the author has indicated that the fundamental 
fluid properties of shear modulus and fluid shear strength can be represented by as few as three terms 
involving pressure and temperature only. 
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Appendix A 

This appendix contains the regression coefficients for the traction data from the two fluids. The 
coefficients for slope are listed under column “m” and the coefficients for the traction coefficient 
are listed under 


Coeffi- 

cient 

Santotrac 50 

Sunoco TDF-88 








P 

u 

m 

P 

P 

m 


eq. (14) 

eq. (15) 

eq. (14) 

eq. (14) 

eq. (15) 

eq. (14) 

U 

44.906 

45.933 

45.890 

22.340 

23.076 

22.476 

P 

1.356 

1.331 

1.369 

1.449 

1.408 

1.424 

Q 

53.491 

53.420 

54.890 • 

49.489 

50.296 

49.076 

K 

3.340 

3.321 

2.753 

2.239 

2.441 

2.262 

Z 


5.786E-03 



.639E-03 


Cl 

0.843E-01 

0.745E-01 

0.450E-02 

0.877E-01 

0.859E-01 

0.452E-02 

C2 

-.685E-02 

-.536E-02 

-.719E-01 

-.515E-02 

-.582E-03 

-.492E-01 

C3 

.434E-01 

-.175E-02 

.114E-02 

.154E-01 

-.480E-02 

.649E-01 

C4 

.421E-01 

-.219E-02 

.109E-01 

-.379E-02 

.581E-03 

-.102E-02 

C5 

-.960E-03 

.204E-01 

-.701E-01 

-.676E-03 

.162E-01 

.202E-01 

C6 

.114E-02 

-.208E-01 

-.817E-01 

-.428E-03 

.211E-03 

-.154E-01 

C7 

-.449E-03 

.109E-01 

.444E-00 

.444E-03 

-.601E-02 

-.559E-00 

C8 

-.120E-02 

-.258E-02 

.805E-00 

-.346E-02 

-.542E-04 

.609E-00 

C9 

-.197E-01 

-.592E-03 

-.755E-01 

-.126E-01 

-.102E-02 

-.124E-02 

CIO 

.299E-02 

-.415E-03 

-.577E-00 

.356E-02 

-.353E-03 

.lOlE-02 

Cll 

.350E-01 

.705E-03 

.124E-02 

.141E-01 

-.136E-03 

.157E-02 

C12 

.664E-02 

.957E-03 

.115E-02 

.764E-02 

.401E-03 

.270E-01 

C13 

.172E-02 

.176E-02 

-.921E-02 

.167E-02 

-.132E-02 

-.639E-00 

C14 

-.317E-02 

.305E-02 

-.759E-01 

.317E-03 

-.848E-03 

.558E-00 

C15 

-.285E-02 

-.588E-02 

.173E-00 

-.851E-04 

-.255E-02 

.245E-01 

C16 

.956E-02 

-.427E-02 

.250E-01 

.532E-02 

.103E-02 

-.196E-01 

C17 


-.162E-01 



-.137E-01 


C18 


-.273E-02 



-.585E-03 


C19 


-.534E-02 



.344E-02 


C20 


-.907E-03 



.143E-02 


C21 


.490E-01 



.277E-02 


C22 


.282E-01 



-.382E-02 


C23 


-.406E-01 



.141E-01 


C24 


-.588E-01 



.701E-02 


C25 


.324E-03 



.177E-02 


C26 


.790E-04 



.104E-03 


C27 


-.138E-04 



-.162E-02 


C28 


-.302E-03 



-.723E-03 


C29 


-.993E-02 



-.850E-04 


C30 


-.112E-01 



.791E-03 


C31 


.212E-01 



.389E-03 


C32 


.202E-01 



-.291E-02 
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